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MODELING THE BROADBAND EMISSION OF FERMI/LAT GRB 090902B 
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ABSTRACT 

GRB 090902B, detected by Fermi Large Array Telescope (Fermi/LAT), shows extend high-energy emission 
(> 100 MeV) up to 10 3 s after the burst, which decays with time in a power-law as t . It has been also 
observed by several follow-up low-energy instruments, including an early optical detection around 5000 s after 
the burst. The optical emission at early time decays faster than t 6 , which has been suspected to originate from 
the reverse shock. We here explore the models that can possibly explain the the broadband afterglow emission 
of GRB 090902B. We find that the reverse shock model for the early optical emission would overpredict the 
radio afterglow flux that is inconsistent with observations. A partially radiative blast wave model, which though 
is able to produce a sufficiently steep decay slope, can not explain the broadband data of GRB 090902B. The 
two-component jet model, which consists of a narrow and bright jet component in the core and a surrounding 
wider and less energetic jet component, is shown to be able to explain the broadband afterglow data, including 
the LAT high-energy data after ~ 50 s and low-energy (radio, optical and X-ray) afterglow data. The early- 
time high-energy emission detected by LAT before ~ 50 s is likely due to internal origin as that of the sub- 
MeV emission. The highest energy (33 GeV) photon of GRB090902B detected at 80 s can be marginally 
accommodated within the forward shock emission under the optimistic condition that electrons are accelerated 
by the Bohm diffusive shock. 

Subject headings: gamma ray: bursts — radiation mechanism: non-thermal 



1. INTRODUCTION 

GRB 090902B, triggered the Fermi Gamma-ray Burst 
Monitor (GBM), is a long, intense burst with a redshift of z = 
1.822 (Cucchiaraet al. 2009) and a fluence of (4.36 ±0.06) x 
10" 4 ergcnT 2 (10 keV-10 GeV) over the duration T 90 ~ 25s 
of the prompt emission. These data give an isotropic energy 
Ei so = (3.63 ±0.05) x 10 54 erg, comparable to the energy in an- 
other bright GRB detected by Fermi, GRB 0809 16C (Abdo et 
al. 2009a). It was also detected by the Fermi LAT with ex- 
tended high energy (> 100 MeV) emission up to 1000 s after 
the trigger, which includes a 33.4 GeV photon at 82s after the 
trigger, the highest energy yet detected from GRBs. The spec- 
trum of this extended emission is consistent with a power law 
with photon index T = -2.1 ±0.1, and its flux (>100 MeV) 
declines as f" 1 - 5±0 1 over the interval from 25 to 1000 s after 
the trigger (Abdo et al. 2009b). 

The optical-infrared afterglow of GRB 090902B has been 
observed by several instruments (Pandey et al. 2010; 
McBreen et al. 2010). The earliest optical detection by 
ROTSE-IIIa occurring at about 5000 s after the burst reveals 
a bright afterglow which decays faster than f" 1 6 (Pandey et 
al. 2010). The temporal decay of the optical afterglow be- 
comes flatter after ~ 12.5 hours with a decaying index of 
-0.90 ± 0.08. The X-ray afterglow, detected by Swift XRT 
(Kennea & Stratta 2009), is consistent with a single power 
law decay with an index of ~ 1.30 ±0.04 from 12.5 hours to 
17 days after the burst. In the radio band, the light curve re- 
mains flat until weeks after the burst and then starts to decline 
(Cenko et al. 2010). 

Kumar & Barniol Duran (2010) modeled the the LAT, X- 
ray and late-time (f > ldays) optical data of GRB090902B 
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with the forward shock synchrotron emission and find that 
one forward shock component can explain these data. How- 
ever, their modeling does not include the early optical data at 
5000 s and the long-term radio afterglow data. Cenko et al. 
(2010) model the late-time optical, X-ray and radio data of 
GRB090902b with one forward shock component. Follow- 
ing Pandey et al. (2010), they attribute the early optical flash 
to the reverse shock emission. They also calculated the ex- 
pected high-energy emission using their best-fit parameters, 
but find that the expected 100 MeV flux at t = 50 s is a factor 
of 7 below the observed value. Given these incompleteness 
and inconsistency, we aim at modeling the broadband data of 
GRB090902B, including the LAT, optical, X-ray and radio 
data. 

One interesting feature of GRB090902B is the fast tempo- 
ral decay of the early optical emission, i.e. the power-law 
decay slope in F v oc t~ a is a > 1 .6. Analogous phenomenon 
has been observed in GRB 990123 and has been interpreted 
as arising from the reverse shock passing through the bary- 
onic ejecta of the GRB outflow, which produces an optical 
flash with a steep decay (typically evolves as f" 2 ) at the early 
time (Meszaros & Rees 1997, 1999; Sari & Piran 1999a, b; 
Panaitescu & Kumar 2004; Nakar & Piran 2004, 2005). The 
reverse shock scenario has been also proposed for the early 
optical emission of GRB 090902B (Pandey et al. 2010; Cenko 
et al. 2010). With the high-energy afterglow data from Fermi 
LAT for GRB 090902B, the forward-reverse shock scenario 
for the all-band data can be further tested. 

There are other scenarios that can possibly explain the early 
fast decay optical and high-energy emission, such as the fully 
or partially radiative blast wave model (Ghisellini et al 20010; 
Wu et al. 2005a). The radiation emitted by electrons accel- 
erated in shocks may lead to energy loss of the shock. The 
adiabatic shock approximation will break down when the en- 
ergy loss becomes important, which occurs when the equipar- 
tition factor for electrons (i.e. e e ) and the radiation efficiency 
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of electrons are high. In the extreme case when almost all the 
shock energy goes into electrons (i.e. e e — 1) and the electrons 
are in the fast-cooling regime, the shock becomes fully radia- 
tive. Although the fully radiative shock is in general unreal- 
istic for GRBs, a partially radiative shock is possible when e e 
is relatively high, especially during the early afterglow phase. 
Due to the decreasing energy remained in the forward shock, 
a faster decay of the afterglow than that predicted by the stan- 
dard scenario is expected (Wu et al. 2005a). Motivated by 
this, we also consider whether this scenario can explain the 
early fast decay of the optical emission in GRB 090902B. 

Numerical simulations show a structured outflow when the 
jet breaks out the collapsing stellar envelope of massive stars 
(Zhang et al. 2003). In the structured jet model, it is usu- 
ally assumed that the energy per unit solid angle depends 
as a power-law or a Gaussian function on the angular dis- 
tance from the axis (e.g. Meszaros , Rees, & Wijers 1998; 
Dai & Gou 2001; Rossi, Lazzati, & Rees 2002; Zhang & 
Meszaros 2002; Kumar & Granot 2003; Granot & Kumar 
2003; Zhang et al. 2004). For the sake of calculation ease, 
the structured jet can be simplified as a two-component jet 
(Ramirez-Ruiz et al. 2002, Berger et al. 2003; Peng et al. 
2005; Huang et al. 2004), which consists of a narrow and 
bright jet component in the core and a surrounding wide and 
less energetic jet component. When observers view along the 
axis, one will see the afterglow emission produced by both the 
bright core and the broad wings surrounding the core. As the 
narrow component has a small opening angle, the jet break in 
the light curve can occur very early and the afterglow emis- 
sion produced by the narrow component will have a fast decay 
after that time. This provides a potential mechanism to pro- 
duce the early-time fast decaying optical emission of GRB 
090902B, while the late radio, optical and X-ray afterglows 
that have normal light curves can arise from the wide compo- 
nent. We will study whether this two-component scenario can 
explain the broadband afterglow data of GRB 090902B. 

We first study the reverse shock scenario in §2, and find 
that the reverse shock scenario for the early optical emission 
predicts a higher radio flux than the observed value. In §3, 
we study whether the partially radiative forward shock can 
explain the broadband data of GRB 090902B, but find a neg- 
ative result. Then we propose a two-component jet model for 
the broadband afterglow data of GRB 090902B in §4. In §5, 
we discuss the origin of early-time high energy (> lOOMeV) 
emission observed during the prompt bursting phase. Finally 
we give our conclusions and discussions in §6. 

2. FORWARD SHOCK - REVERSE SHOCK (FS-RS) MODEL 

When the ultrarelativistic cold baryon-dominated GRB 
ejecta encounters the cold ISM, a reverse shock that propa- 
gates back into the ejecta and a forward shock that propagates 
into the ISM will form, and as a result both the ejecta matter 
and the ISM matter are heated up. An optical flash is expected 
to be produced by the reverse shock synchrotron emission, 
which should decay quickly due to that the shocked ejecta is 
expanding adiabatically. The early optical flash from GRB 
990123 with a 9th magnitude is believed to be such a good 
case (Sari & Piran 1999a; Meszaros & Rees 1999). It re- 
mains unclear why such optical flashes are lacked in general 
in GRBs as ground-based robotic telescope observations (Yost 
et al. 2007) and early Swift UV/Optical Telescope (UVOT) 
observations (Roming et al. 2006) yield non-detections down 
to a much lower limit. 

A bright optical flash is detected by ROTSE from GRB 



090902B with a magnitude of m R = 16.4 ±0.5 at t = 5320 
s after the burst. The optical non-detections at 2 x 10 4 s by 
ROTSE imply that the optical emission decays faster than 
r L6 . Pandey et al. (2010) and Cenko et al. (2010) argue 
that this rapid decay slope is suggestive of the reverse shock 
origin. With the rich multiwavelength observation data avail- 
able for this burst, including the high-energy LAT data, the 
radio data, late-time optical and X-ray data, we aim at testing 
this reverse-shock origin possibility. 

2. 1. The forward shock emission 

The late-time optical and X-ray afterglow emission after 
half a day shows a normal decay, with a decay slope of 
a = -0.90 ± 0.08 and a x ~ -1 .30 ± 0.04 respectively. They 
are broadly consistent with the synchrotron afterglow emis- 
sion produced by a forward shock expanding into a constant- 
density medium with an electron index of p ~ 2.2 if the cool- 
ing frequency in the synchrotron spectrum is between the op- 
tical and X-ray frequencies (Kumar & Barniol Duran 2010; 
Cenko et al. 2010). Following Sari et al. (1998) and Wijers & 
Galama (1999), we get the radius of forward shock at time T 
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where E is the isotropic-equivalent kinetic energy of the 
burst and n is the number density of the circum-burst ISM. 
Throughout the paper, and we use c.g.s units and denote by 
Q x the value of the quantity Q in units of 10*. Adopting 
R = 4--f 2 cT, the bulk Lorentz factor of the shocked matter can 
be described by 
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Then we can obtain the characteristic frequencies of the after- 
glow synchrotron emission 

X /2 4?- 5 (l+z) V V /2 Hz (3) 



■ 1.05 x 10 18 f 2 (p)e 2 e 



where f(p) = 6(p-2)/(p- 1), e e f and esf are the equipartition 
factor for the energy in electrons and magnetic field in the 
shock, and 

v c t = 1.13 x io 24 £ 5 " 4 1/ V^B 3 f (_ 2 5 [i+y(^ cf )]- 2 (i+z)- 1 / 2 r " 1/2 Hz 

(4) 

where F(^ c f) is Compton parameter for electrons that pro- 
duce the synchrotron photons at v c f. Hereafter the super- 
scripts/subscripts 'f and 'r' are used to represent the quanti- 
ties of the forward shock and reverse shock respectively. The 
peak flux density is 

^ f max = 4.08 x \Q- 2 \% 2 E 5 ^ 5 {\+z)D-^xgcm-h- x (5) 

Then we get the flux density at a fixed frequency 



(v/v m fT (p - 1)/2 oc T~ 3 ^ 4 

(^mfr ,V W^mf) 1/2 

cxr- ( >- 2) / 4 [i+y( J /)]- 1 



V < fmf < ^cf, 
V mf < V < Ud, 

fmf < fcf < V. 



(6) 

where Y(v) is the Compton parameter for the electrons whose 
synchrotron frequency is v. Due to Klein-Nishina scattering 
effect, Y(v) is not a constant for the electrons that produce 
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high energy photons (Wang et al. 2010) and in fact is depen- 
dent of v. 

Now we confront the above theory with the observed data to 
constrain the unknown parameters. From Abdo et al. (2009), 
Pandey et al. (2010), we take the observed flux for LAT data at 
173s, R-band optical data at 1 .30 x 10 5 s, X-ray data at 1 .09 x 
10 5 s and radio data at 4.67 x 10 5 s as below: 

F/ AT (173s) ~ 3.74 x lO-^ergcnrV 1 

(l/LAT > V cf > V mf )) 



F|(1.09 x 10 5 s) ~ 1.82 x lO-^ergcnrV 1 

(^X > ^cf > ^mf) 

KF o f pt (1.30 x 10 5 s) ~ 9.33 x lO-^ergcnrV 1 

(^cf > ^opt > fmf) 

(4.67 x 10 5 s) ~ 6.83 x lO-^ergcnrV 1 

(fmf > t'radio > V a f) 



F f 

radio 



where n ~ 0.85 is the extinction coefficient of the host galaxy 
in R band (Cenko et al. 2010). The spectral regime of 
each frequency is given in the parentheses. Although z^lat 
and vx belong to the same spectral regime, the X-ray flux 
and high-energy flux data can give two independent con- 
straints due to two different Y parameters. For electrons 
that produce the X-ray synchrotron afterglow, the inverse- 
Compton scattering loss is in the Thompson regime and Y ~ 

74.3e e 2 £ ieB f 5J E 5 1 4 /18 ny 3 18 (l+z) 1 / 18 r - 1/18 in the slow-cooling 
phase, while for high-energy gamma-ray photons, Y parame- 
ter is small due to the KN suppression effect on the scattering 
cross section, as we will show later. 
So we obtain the constraints on the shock parameters 
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12.9[l+T(100MeV)] 

51.6 
- 13.5 
0.952 



(8) 



where y(100MeV) is the Compton parameter for the electrons 
whose synchrotron frequency is hv = lOOMeV at t = 173s. In 
the calculation, a typical value of p = 2.2 for forward shock 
has been used. Solving the above equations, we obtain 

' £54 ~2.0[l+T(100MeV)] 3 / 4 
«_ 3 ~ 0.40[1 +T(100MeV)] 1 / 4 
e ef _i ~4.2[l+T(100MeV)] 1 / 4 ( ) 

k e Bf -5 ~9.2[l+T(100MeV)]- 7 / 4 . 

For the above set of parameters, we obtain T(lOOMeV) ~ 

Qmf{p)-'lh^y B f 5 Ell\%XI\l + zyi" ~ 0.06 at t = 
173s (Wang et al. 2010). As T(lOOMeV) < 1, we can just 
omit T(lOOMeV) in these expressions. 



2.2. The reverse shock emission 

As the reverse shock propagates into the fireball ejecta, the 
number of shocked electrons increases, which leads to an ini- 
tial rise of the reverse shock flux. Once it has passed through 
the shell, no new electrons will be shocked and the flux will 
decrease due to adiabatic expansion of the radiating gas. The 
flux of reverse shock peaks at T; = max(7d ec , 7c r0S s) (hereafter, 
we use the subscript T to represent quantities at the peak 
time). Here 7d ec = (//r/ 8//3 c)(l +z) is the shell deceleration 
time, where I = (E /nnipC 2 ) 1 / 3 is the Sedov length and 77 is 
the initial Lorentz factor of the ejecta. r cross = (A/c)(l +z) = 



T90, is the characteristic timescale within which the reverse 
shock crosses the shell, where A is the thickness of the shell 
(Sari & Piran 1995, 1999b). When the shell is thick so that 
Tcross > Tdec, the reverse shock becomes relativistic early on, 
and we call it "thick shell" case, otherwise it belongs to the 
"thin shell" case. Evolution of bulk Lorentz factor, gas pres- 
sure and density of shocked electrons after crossing time fol- 
lows 7 oc r- 7 / 16 , p oc T" 13 / 12 and n oc r -13 / 16 (Sari & Piran 
1999a, b; Meszaros & Rees 1999; Kobayashi 2000). 

The bulk Lorentz factor of the shocked shell at the peak 
time is 7 = 1838(1 +z) 3 / 8 £ 5 1 4 /8 n: 3 /8 r i " 3/8 . Then we can get the 



(7) Lorentz factor of the reverse shock at the peak time 

- 1 ( n _,_ 7i 
7 = x — + — 

2 \7< V 
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The minimum random Lorentz factor of electrons is given by 
p — 2m„ 

7mr = £er T~ (7~1) 

p-1 m e 



(11) 
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where e er is the fraction of thermal electrons carried by mag- 
netic field in reverse shock, b = Szil (namely, 7 — 1 = b^) 

and i=T/Ti. We use this approximation for 7 - 1 so that we 
can treat the calculation analytically. Usually, the value of b 
is taken to be 1 (Waxman & Draine 2000; Kobayashi et al. 
2007) or i (Wang et al. 2005). A mildly relativistic reverse 
shock approximation with 7 - 1 ~ 1 is also sometimes used 
(see e.g. Sari & Piran 1999a). In this paper, we take b = \ in 
the calculation, as r\ > 7,. For a typical value of p = 2.5 for 
the reverse shock emission, f(p) = 2. 

We get the two characteristic frequencies for the reverse 
shock synchrotron emission 

v m , = 9.24 x 10 1 W(/^3 2 4-i4f-i"-3 2 (l +z)" 1 r 73 / 48 Hz 

(12) 

-1/2,1/16, 



and 



i.i2xio 18 £ 5 - 4 1/2 «:U: 3/2 ' 



(13) 

The peak flux of the reverse shock synchrotron emission is 



1 v, max 
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=9.75 X lO" 22 ^. 
(1 +z) 7 / 4 DL 2 28 r i - 3 / 4 r 47 / 48 ergcm- 2 s- 1 . 
Then we get the flux density at a certain frequency 



(14) 



F t =F x 

v ZAmax 



(^/Vmrr ( "- 1)/2 OC r-( 7 > +21 )/ 96 Vrt < V < U a 



(15) 

If v > max(i' mr ,;/ cr ), the flux drops exponentially with time 
because there are no new injected electrons whose typical syn- 
chrotron frequency lies in this regime. 

If we attribute the early optical emission at t ~ 5000 s to the 
reverse shock synchrotron emission, we would have 



KF o r pt (5000s) = 0.91mJy 



(16) 



where n ~ 0.85 is the extinction coefficient of the host galaxy 
in R band (Cenko et al. 2010). The reverse shock emission 
can produce a radio flare as has been seen in GRB 990123 
(Kulkarni et al. 1999). Although the radio emission in GRB 
090902B drops after the first detection, it shows a flattening 
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after the second detection at 4.67 x 10 5 s, so the forward shock 
emission should be dominated at this time. Thus, we expect 
that 



^(4.67 x 10 5 s) < ^ (4.67 x 10 5 s) 

= 3.41 x 10" 28 ergcm" 2 s _1 . 



(17) 



At 5000 s, from equation (12) and (13) and using the con- 
strained value of E and n from the forward shock emis- 
sion, we get v mx ~ 6.39 x 10 8 77 2 e 2 T _ 1 e B (_ 1 Hz and v C! ~ 3.29 x 

lO'^g^Hz, so the optical frequency ^ opt (4.7 x 10 14 Hz) is 
expected to be in the regime v mr < v opt < v a and the flux den- 
sity should decline with a temporal index of ~ -2. If the radio 
frequency z^ ra dio (8.5 x 10 9 Hz) also lies in the same regime (i.e. 
^mr < ^radio < ^cr) at 4.67 x 10 5 s, the radio flux then relates 
with the optical flux as 



1 radio 



5000s 



8.5 x 10 9 Hz J 



Combining Eqs.(17) and (18), we obtain p < 1.5. Such 
a small p is inconsistent with the observed decay slope 
of the late-time optical and X-ray emission. As v ml oc 
f~ 73 / 48 , we note that at 4.67 x 10 5 s, v w > f rad io only if 773 > 
1 12c" 1 ..! e B J.^i • Since such an initial Lorentz factor is too large, 
this spectral regime is unlikely (see c.f. Ioka 2010). 

Radio flux from the reverse shock may be also affected by 
the synchrotron self absorption (SSA) of the radiating elec- 
trons. The SSA frequency in the slow cooling case is (Wu et 
al. 2005a) 



a 3 / 5 tW,Z/ar < ^mr, 

a 2/(p+4) v mr7 v mr < v w , 



(19) 



where 

_ c (p-l)eT, 



a = 



3777 ? 3 6 J B54n:re; r 5 - 1 e B ,-i(l +zfTr 2 f 



(20) 

co ~ 15 is nearly a constant and £ = E ^^i is the column 
density of electrons heated by the reverse shock. The SSA 
effect leads to a different light curve in the radio band from 
that in the optical band if f ra dio lies in one of the following 
regimes: z/ radio < v mx < v- dI ,v ml < f rad io < ^ ar , and v mdio < v w < 
v ml . Among these three cases, the second one is the most 
likely case for GRB 090902B since we expect v mr < f ra dio 
for normal parameter values. So we just consider this case 
here and leave discussions on other cases in the Appendix. 
We find v mx < v m &io < v. M requires £ , 5 4n_3e er _ie B r-i > 10 4 at 
4.67 x 10 5 s according to equations (3), (19) and (20). If we 
substitute £54 = 2.0 and n_3 = 0.40, which are derived from 
the constraints by the forward shock emission, into this in- 
equality, we obtain e er -ieBr,-i ^ 1-25 x 10 4 , which is unlikely. 
Therefore we conclude that SSA effect can not solve the prob- 
lem that the reverse shock scenario for the optical emission 
overpredicts the radio emission. 

3. PARTIALLY RADIATIVE MODEL 

It is usually assumed that the blast wave that produces 
the afterglow emission is nearly adiabatic (Meszaros & Rees 
1997; Sari et al. 1998), i.e. the total kinetic energy of the 
relativistic shock is a constant. However, the energy loss of 
the blast wave could be significant under some circumstance. 



The radiation efficiency of the blast wave is given by (Wu et 
al. 2005a) 



/ \ (p-2)/2 



(21) 



As we can see from the above equation, the energy loss is 
especially important in the "fast-cooling" case (i.e. v m > v c ) 
with a large e e . In this case, the decreasing blast wave en- 
ergy at early times will result in a faster decay of the after- 
glow emission than the adiabatic case. At late time, as — 
decreases with time, the radiation efficiency drops and the de- 
cay slope changes to the adiabatic case. We examine whether 
this scenario can explain the fast decay in the early optical and 
high-energy LAT emission. 

According to Huang et al. (1999) and Wu et al. (2005a), the 
isotropic-equivalent energy E of the blast wave evolve with 
time as 

, j x -3e/(4-e) 

E ' E i^) ■ <22) 

The quantities describing the synchrotron spectrum in such 
a semi-radiative shock are similar to equations (3)-(5), ex- 
cept that the constant E in these equations should be replaced 
by a time-dependent E as described by Eq. (22). So the 
synchrotron emission flux decays as F v oc 7 , - 3 0>- 1 +<0/( 4 -<0 f or 
v m { < v < v c f. To explain the decay slope a > 1 .6 of the early 
optical emission in GRB 090902B with this model, we need 
e > 0.6. Thus we need 



£ef>0.6. 



(23) 



To explain the late optical, X-ray and radio observations, we 
also need: Kf f p t( 1()5s ) - 0-OlmJy, F^IO^) ^ 0.03mJy, and 
Fx(10 5 s) — 0.2/xJy (radio observations starts at about 1.3 x 
10 5 s, we extrapolate it to 10 5 s). At 10 5 s, the two character- 
istic frequencies in the synchrotron spectrum are v m { ~ 2.0 x 

lO^fe^Hz and v ci ~ 1.1 x lO 1 ^ 7 V^^Hz 
for e e = 0.6 and p = 2.2. As v opt , vx, Vm&io lie in three different 
frequency regimes, we have three independent constraints on 
the shock parameters. Finally we obtain 



n_ 3 ~ 0.53,e Bf .- 5 ^ 0.39,£ 54 (10 5 s) ~ 8.0 



(24) 



For these parameters, the deceleration time of the blast wave 
is Tdec — 120s, and we can obtain the blast wave kinetic en- 
ergy at the deceleration time Tdec according to Eq.(??), i.e. 
£54(7dec) — 50. This energy is extraordinary large. With 
such a high isotropic energy, we expect the flux density in 
LAT band at the deceleration time to be FL AT (7de C ) ~ 1/xJy, 
which is one order of magnitude higher than the observed flux 
(~ 0.1 /xJy). Therefore, we conclude that this model can not 
explain the broadband data of GRB 090902B. 

4. TWO-COMPONENT JET MODEL 

Jets from GRBs may have complex structure. For the sake 
of calculation ease, the structured jet can be simplified as a 
two-component jet. It assumes that the jet consists of two 
components: a narrow component with a relatively small half- 
opening angle (9 N ) and a large isotropic-equivalent energy in 
the center, and a wide component with a larger half-opening 
angle (fiw) and a smaller isotropic-equivalent energy (here- 
after, we use the superscripts/subscripts 'N' and 'W' repre- 
sents the quantities of the narrow component and the wide 
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component respectively). The energy distribution of the two- 
component jet can be parameterized as: 

p ia\-l Eiso.N <9 <9 N , fr)t .. 

tisoW- 1 E . soW e N <e<e w {z:>) 

The two components of the jet encounter the ISM and gener- 
ate forward shock and reverse shock respectively. The ob- 
served afterglow emission is the superposition of the two 
component. Due to a higher kinetic energy, the contribution 
by the narrow component to the afterglow emission should 
dominate at the early time. As the bulk Lorentz factor 7^ of 
the narrow component decreases, the light curves breaks to a 
steeper one once 7^ < 1/9^. After the break, the light curve 
of the narrow component declines much faster than that of the 
wide one, and the contribution by the wide component could 
dominate at late times. For the case of GRB 090902B, we 
can attribute the early optical emission around 5000s to the 
forward shock emission of the narrow component after the 
jet break, since the optical emission decays very fast. The 
extended high-energy emission detected by LAT can be at- 
tributed to the afterglow emission of the narrow component 
before the jet break and the late radio, optical and X-ray emis- 
sion can be attributed to the afterglow emission of the wide 
component. 

Considering the inverse Compton loss of electrons, the syn- 
chrotron flux density at frequency v after the jet break is 
(Rhoads 1999; Sari, Piran & Halpern 1999) 

' (*V« 1/3 oc r-'/ 3 

(v < v m < v c ), 



(v/v m T p/2 (vc/v m ) 1/2 oc T-"[l+Y(u)T l 

Om < V c < V). 

(26) 

The jet break of the narrow component should occur after the 
last LAT detection (around 1000 s) and before the first ROTSE 
optical data at about 5000 s. Assuming a jet break time at, e.g. 
T Break ~ 4000s ' the broadband data of GRB 090902B imply 

' F^ T (173s) ~ 3.74 x lO^'ergcrrrV 1 

kF o *(4000s) ^ 9.06 x lO-^ergcnrV 1 

< KF f t (1.3 x 10 5 s) ~ 9.33 x lO^ergcrrrV 1 (27) 

Ff (1.09 x 10 5 s) ~ 1.82 x lO^ergcm-V 1 
> F ZJ 4 - 61 x lo5f O - 6 - 83 x lO-^ergcnrV 1 

By comparing the above flux data with the theory prediction 
of the two-component jet model, we get the constraints 

' 4i-^m-5 E UN - 12.9[l+7(100MeV)] 

6/5 4/5 13)10 1/2 
e ef,-l t Bf.-5 c 54.N "-3 — 
, 6/5 4/5 £,13/10 1/2 , ,~ S s 

' £ ef,-l £ Bf.-5 £ 54.w"-3 - 5L6 ( 28 > 

8/15 1/2 „ 1/2 ,„ - 

e ef _ieBf-5 £ 54>W«_ 3 — 13.5 

-2/3 1/3 5/6 w l/2_n 952 
I. fc ef -1 fc Bf -5 C 54,W -3 — U - 5 ' JZ - 

Solving this set of equations, we obtain 

' £ 5 4.n ~2.7[l+F(100MeV)] 3 / 4 

F 54 *w ^ 0.65[1 +F(100MeV)] 3 / 4 
< n_ 3 -O^tl+^lOOMeV)] 1 / 4 (29) 
e ef _i ~2.9[l+T(100MeV)] 1 / 4 

k e Bf .-5 ^ 130[l+F(100MeV)]- 7 / 4 . 



100Mev Emission 




Time since burst (s) 

Figure 1. Fit of the light curve of the high-energy (> lOOMeV) gamma-ray 
emission from GRB 090902B with the two-component jet model. The dashed 
line and dash-dotted line represent the forward shock synchrotron emission 
of the narrow component and wide component respectively (note that the two 
lines overlap during the rising phase), while the solid line is sum of the con- 
tributions by two jet components. The green, red and blue lines correspond to 
the initial Lorentz factor of 3000, 1800 and 1300 respectively. The black dot- 
ted lines represent the high-latitude emission of the last prompt LAT pulse, 
calculated with Eq.(33), assuming pulse variability time of AT =0.05s, Is, 
5s from left to right respectively. Data are taken from Figure 2 of Abdo et al. 
(2009b). The fitting parameters are shown in Table 1. 



Optical Emission 




Figure 2. Fit of the light curve of the optical emission from GRB 090902B 
with the two-component jet model. The thick dashed line and dash-dotted 
line represent the forward shock synchrotron emission of the narrow com- 
ponent and wide component respectively, while the solid line is sum of the 
contribution by two jet components. The thin dashed line and dash-dotted 
line represent the reverse shock synchrotron emission from the narrow com- 
ponent and wide component respectively. Data are taken from Figure 6 of 
Cenko et al. (2010) and Figure 11 of McBreen et al. (2010). The parameters 
used in the calculation of the reverse shock emission is shown in Table 2. 



For the above parameter values, we find that F(100MeV) <C 1 
(Wang et al. 2010). 

Figures 1-4 show the fit result of the LAT, R band, X- 
ray and radio data of GRB 090902B, respectively, and the 
fitting parameters are given in Table 1. In the LAT en- 
ergy band (Fig. 1), we plot the forward shock emission 
for three different values of the initial Lorentz factor, i.e. 
T] = 3000, 1800 and 1300 respectively. In the case of r] = 3000, 
the shell crossing time r cross = Tgo = 25s is larger than the 
shell deceleration time 7d ec , so the reverse shock belongs to 
the thick shell case. In this case, the reverse shock will 
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X-ray Emission 




Figure 3. Fit of the light curve of the X-ray emission from GRB 090902B 
with the two-component jet model. All the lines represent the same as those 
in Figure 2. Data are taken from Figure 6 of Cenko et al. (2010). 



Radio Emission 




10-" I , , , , i , 

10 4 10 5 10 s 10' 

Time Since burst {s} 

Figure 4. Fit of the light curve of the radio emission from GRB 090902B 
with the two-component jet model. All the lines represent the same as those 
in Figure 2. Data are taken from Figure 6 of Cenko et al. (2010). 

transit from the Newtonian phase to the relativistic phase 
at time t N = = 0.38£ 5 'f nl3 /2 (r//3000)- 4 (l +z) 3 / 2 s. At 

time tff < t < r cl0SS , the Lorentz factor of the forward shock 
evolves as T oc f -1 / 4 , and the radius of the shell evolves as 
R oc t x l 1 . Thus, the characteristic frequencies and the peak 
flux evolve as v m f oc t~ l , oc f" 1 and Fl max oc t . As a result, 
the flux of high energy emission evolves as Flat oc f H>+3)/2 in 
the regime v m f < ^lat < Vcf> an d evolves as Flat oc fH>+2)/2 
in the regime ^lat > v c f. For the inferred parameters of 
GRB090902B, </ L at > ftf, so Flat oc fH>+2)/2 K f -o.i before 
Tci-oss (for p = 2.2). We can see that the rj = 3000 case (the 
green solid line) can not account for the flux of the LAT peak 
at 10-20 s, although its light curve peaks at the right time. 
When r\ < 1800, the reverse shock belongs to the thin shell 
case and light curve peaks after Tgo- In the thin shell case, the 
light curve rises initially as f 2 before the peak. For ij = 1800, 
the forward shock emission can account for the LAT data after 
~ 50 s (the red solid line), but falls below the observed high- 
energy flux before ~ 50 s. When r/ is smaller, the deceleration 
time is longer, as seen by the the case of 77 = 1300 (the blue 
lines). 



In the optical band (Fig. 2), the forward shock emission of 
the narrow component reproduces the early optical flash at 
~ 5000 s. The isotropic energy in the narrow jet compo- 
nent is £N,iso = 2.7 x 10 54 erg and its half opening-angle is 
9n = 0.36°. The afterglow emission of the wide component 
can also reproduce the late-time (t > 1 days) optical data as 
well as the X-ray and radio data, if there is a jet break at ~ 6 
days (see Figs 2-4). The discrepancy between the model light 
curve and the first radio detection at ~ 10 5 s could be due to 
strong interstellar scintillation of radio emission at early times 
(Cenko et al. 2010). In the fitting, the half opening-angle of 
the wide component is 9^ = 2.8° and the isotropic energy is 
£\v,iso = 0.65 x 10 54 erg. The circum-burst density in the fitting 
is n = 0.28 x 10" 3 cm" 3 , which is at the low end of the ISM 
densityQ Our fitting parameters of the wide component agree 
well with that of Cenko et al. (2010), who modeled only the 
late-time optical, X-ray and radio afterglow emission. With 
the extra contribution by a narrow jet, our two-component jet 
model is able to fit also the early optical data and the LAT 
high-energy data after ~ 50 s. Due to a relatively low e e and 
a low circum-burst density n that lead to a low radiative effi- 
ciency, the radiative energy loss of the shock is unimportant, 
so the adiabatic shock approximation used in the calculation 
is valid. 

Ryde et al. (2010) argue that a thermal photosphere compo- 
nent is seen in the prompt burst emission of GRB 090902B. 
If this is true, the composition of the jet should be baryon- 
dominated. In this case, the reverse shock emission produced 
by the two components should also be present. We calculate 
the reverse shock emission produced by the two component 
jet, which are also shown in Figures l—^. We find that the 
reverse shock emission is typically subdominant and that the 
parameters of the reverse shock are not well constrained. The 
parameter values given in Table 2 are just one illustration. 

5. ORIGIN OF THE EARLY-TIME HIGH-ENERGY (> 100MEV) 
EMISSION 

We have shown that the forward shock emission of the nar- 
row jet component can account for the long-lived high-energy 
emission after ~ 50 s, as shown in Fig. 1 , but it can not account 
for the early-time high-energy data. The early LAT data is 
inconsistent with a forward shock origin in the following as- 
pects: l)that the high-energy emission of GRB090902B starts 
to decay before is unexpected because before 790 the shell 
is putting energy into the forward blast wave continuously to 
keep it from being decelerated quickly; 2) the LAT peak at 
10-20 s is very sharp and the decay slope immediately after 
the peak is too steep while the forward shock emission model 
predicts a smooth and round peak; 3) the temporal variability 
of the LAT emission during the prompt bursting phase is cor- 
related to some extent with the low-energy emission detected 
by GBM, while the the forward shock emission model pre- 
dicts a smooth light curve. Therefore, there must be an extra 
component that produces the high-energy emission before ~ 
50 s in GRB 090902B. The same situation is seen in the mod- 
eling of GRB090510 (He et al. 2010), in which the forward 
shock emission can not account for the high-energy emission 
before ~ 3 s. 

3 Low density circum-burst environment has also been found for 
other bright Fermi/LAT GRBs such as GRB080916C, GRB090323 and 
GRB090328 (Cenko et al. 2010). The reason why such bright Fermi/LAT 
GRBs have preferentially low circum-burst density is unknown. See Cenko 
et al. (2010) for discussions on the possible explanations for such a low den- 
sity in these bursts. 
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We first check whether the reverse shock can produce the 
peak at 10-20 s in the LAT light curve. Due to that the flux 
of reverse shock synchrotron emission rises rapidly with time 
before reaching its peak and that the reverse shock operates 
only once, the light curve of the reverse shock self inverse- 
Compton emission is expected to form a peak at the crossing 
time (Wang et al. 2001a,b; Kobayashi et al. 2007). However, 
due to the low density inferred for GRB 090902B from the 
broadband afterglow data, the scattering optical depth of elec- 
trons in the ejecta is low, so the IC flux is found to be weak. 
For the parameter values given by Eq.(29), the scattering opti- 
cal depth in the reverse shock ejecta of the narrow component 
at the crossing time is 

a rNl,N 01 1A _7 -l,-,l/2 1/2., , .1/2^-1/2 
Tr ' N = ~4^~ = % £ 54,N"-3 (! +Z) ' T i (3Q) 

« 6.3 x 10 -i V 

where A^ N = E^/r]m p c 2 is the number density of electrons in 
the narrow component jet. According to Eq.(14), the peak 
flux of the reverse shock synchrotron emission is F^" N w 

7.3773 1 eg( 2 _ 1 Jy, then we get the peak flux of the SSC emission 
from the reverse shock 

O = ^S~ - 46 % 2 4(>Jy- (3D 

Since this peak flux is lower than the observed high-energy 
flux, which is ~ 3/iJy (Abdo et al. 2009b), and the peak fre- 
quency of the SSC emission locates at low energy, the SSC 
emission from the reverse shock can not account for the LAT 
peak. Using the parameters constrained by the low-energy 
broadband afterglow data, we find that the synchrotron flux of 
the reverse shock at the LAT band is of the order of 10" 2 /zJy, 
which is also much lower than the observed value. 

Therefore, the early LAT emission has to be attributed to 
some internal dissipation process that occurs at radius much 
smaller than the deceleration radius of the external shock. We 
note that there are pulses around 7.0s, 7.8s, 8.3s, 9.6s and 16s 
in LAT band, most of which have corresponding GBM pulses. 
The coincidence supports the viewpoint that the early-time 
LAT emission is due to internal dissipation within the outflow, 
such as internal shocks. Ryde et al. (2010) show that the time- 
resolved spectra of the prompt emission of GRB090902B 
from KeV to GeV can be decomposed into two components, 
one is the thermal multi-color blackbody component peaking 
at sub-MeV and another is a power-law component extending 
into GeV. Within this picture, the high-energy emission is at- 
tributed to the non-thermal emission (synchrotron emission, 
SSC or Comptonization of the thermal photons) produced by 
electrons in internal shocks (e.g. Ryde et al. 2010; Toma 
et al. 2009; Pe'er et al. 2010). Hadronic scenarios have 
also been proposed, including the proton synchrotron radia- 
tion (Razzaque et al. 2010) or photohadronic cascade radi- 
ation (Asano et al. 2009), which usually need a very large 
energy budget due to the low radiation efficiency of protons 
(Wang et al. 2009). 

The highest-energy (33 GeV) photon is detected after the 
prompt burst, at 80 s after the trigger. In its local redshift 
(z = 1 -822) frame, this energy is 94 GeV. Whether such a high- 
energy photon can originate from the synchrotron radiation of 
the forward shock is an interesting issue (Abdo et al. 2009b; 
Barniol Duran & Kumar 2010; Piran & Nakar 2010), as the 
maximum energy of electrons is limited by shock accelera- 
tion. At 80 s, the bulk Lorentz factor is 7 ~ 700 for the pa- 



rameter values constrained by the broadband afterglow data 
of GRB090902B. So the maximum synchrotron photon en- 
ergy is 

e max ~ 40 K - 1 ( 7 /700)GeV (32) 

according to Eq.(12) in Wang et al. (2009), where n > 1 is 
a parameter describing the efficiency of shock acceleration 
with k = 1 corresponding to the fastest shock acceleration — 
the Bohm diffusive shock acceleration in which the scatter- 
ing mean free path equals to the particle gyroradius. This 
means that the highest energy photons of GRB090902B can 
be marginally accommodated by the forward shock emission 
under the optimistic condition that particles are accelerated 
by the Bohm diffusive shock. Another possibility is that this 
highest energy photon belongs to the prompt component, al- 
though the detailed radiation mechanism is unknown. In the 
latter scenario, the deceleration time of the forward shock 
should be later than 80 s, such as the r\ = 1300 case (the blue 
lines) shown in Figure 1 . 

Now we study the origin of the steep decay immediately 
after the LAT light curve peak at t = 10-20 s. A possible 
scenario for the steep decay of the high energy photons is the 
high-latitude emission of the jet at the end of the prompt emis- 
sion phase. Because photons from high latitude regions with 
respect to the line of sight will arrive later than that from low 
latitude region due to the curved front surface of the jet, one 
observes a fast decreasing emission from the high-latitude re- 
gion rather than an abrupt stop of the emission. The flux den- 
sity of the high-latitude emission evolves with time as (Kumar 
& Panaitescu 2000; Wu et al. 2006; Toma et al. 2009) 

F V {T) = F V (T C ) [ (33) 

where AT is the duration of the last prompt pulse detected 
by LAT, T c is the LAT peak time and (3 is the spectral index 
of high-energy photons. According to Figure 1 of Abdo et al. 
(2009b), there is a LAT pulse around 16s. The variability time 
in the energy band above 100 MeV is not well-determined, 
although the sub-MeV emission detected by GBM shows 
variability timescale of ~ 0.05 s. In some internal dissipa- 
tion models, such as the residual internal shock model in 
which high-energy emissions arise from much larger radii (Li 
2010), the variability time in high-energy emission could be 
longer than that of the low-energy emission. Thus, we take 
/? = 1.1, T c = 16s and assume three different values for AT 
(i.e. Ar =0.05s, Is and 5s respectively) to test the high- 
latitude emission model. We calculate the high-latitude emis- 
sion according to Eq.(33), which is shown by the black dotted 
lines in Fig.l. As we can see, only with a very large AT (i.e. 
AT > 5s) can the peak-time data be fitted by this model. As 
such a long variability is inconsistent with the LAT data of 
GRB090902B, the steep decay may just reflect the tail of the 
on-axis prompt emission. We note that a similar conclusion 
has been reached for the early steeply decaying high-energy 
emission in GRB090510 (He et al. 2010). 

6. CONCLUSIONS AND DISCUSSIONS 

With the broadband data from radio frequencies up to the 
LAT high-energy (> lOOMeV) band available, GRB 090902B 
is a good case to examine the origin of the high-energy emis- 
sion in GRBs detected by Fermi/LAT. This burst has a bright 
optical flash detected by ROTSE at ~ 5000 s, which has been 
suspected to arise from the reverse shock (Pandey et al. 2009; 
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Cenko et al. 2010). We first try to fit the broadband after- 
glow observations of GRB 090902B with the FS-RS model, 
with the reverse shock emission explaining the optical flash 
detected at 5000 s. We find that the optical and radio ob- 
servations cannot be explained simultaneously by this model 
because a bright reverse shock optical emission will yield a 
much higher radio flux than observed. The self-absorption 
frequency of the reverse shock emission is found to be below 
the radio band under the constraints by the late-time forward 
shock optical and X-ray emission, so the self-absorption ef- 
fect can not suppress the reverse shock radio emission. 

Considering that the partially radiative blast wave scenario 
can induce a fast decay of the afterglow emission, we further 
test the partially radiative blast wave scenario, but find that 
it can not explain the broadband data of GRB 090902B ei- 
ther. This is mainly because that late-time optical and X-ray 
flux constraints in combination with the large energy loss in 
this scenario lead to a huge initial kinetic energy in the blast 
wave which overpredicts the early-time high-energy gamma- 
ray emission. 

Then we propose that the two-component jet model, which 
consists of a narrow and bright jet component in the core and 
a surrounding wide and less-energetic jet component, is able 
to account for the broadband observations of GRB 090902B. 
The long-lived high-energy emission and early time optical 
emission can be attributed to the forward shock synchrotron 
emissions of the narrow component. The first optical detec- 
tion at 5000 s should be later than the jet break time of the nar- 
row component so that a fast decay optical afterglow is seen. 
From this, we derive that the half-opening angle of the narrow 
component is 9^ ~ 0.36°. On the other hand, the late-time op- 
tical, radio and X-ray afterglows can be attributed to the wide 
component. To model the radio, X-ray and late-time optical 
light curve, ajet break at ~ 5-10 days is needed. From this, we 
derive the initial half opening angle of the wide component, 
9w — 2.8°. The probability of observing within the central 
bright core is only ~ 10~ 2 . This is consistent with the rare de- 
tection of bursts with energy as large as that of GRB090902B. 

When observers view along the axis, the burst would be 
bright and one will see the afterglow emission produced by 
both the bright core and the broad wings surrounding the 
core. However, when observers view the burst off-axis, one 
would miss the bright core and will see the afterglow emis- 
sion produced dominantly by the quasi-uniform wing. As 



a consequence, we expect that the two-component jet struc- 
ture can be discerned more easily in brighter bursts, such as 
GRB 090902B, through afterglow observations. We note that 
the two-component jet model is favored in another strong, 
long burst - GRB 08031 9B (Racusin et al. 2008). A two- 
component model has also been invoked to explained the 
broadband data of a short GRB 090510 (Corsi et al. 2010), 
with the wide component explaining a mild excess in the op- 
tical band at late times 0. Extended high-energy emission is 
especially useful to diagnose the narrow component, since at 
such early time, high-energy emission is predominantly pro- 
duced by the narrow core component. 

The early LAT emission of GRB090902B before - 50 s af- 
ter the trigger can not be explained by the external shock emis- 
sion and should be due to an internal origin. This is consistent 
with the multiple-pulse structure of the high-energy emission 
and its temporal correlation with the sub-MeV emission seen 
during the prompt bursting phase. Modeling of the broadband 
data of GRB090510 also show that the external shock emis- 
sion cannot account for the high-energy emission before ~ 3 
s after the trigger (He et al. 2010). The fact that one single 
Band function component can fit the spectrum of the prompt 
emission from 10 keV to GeV in GRB080916C also supports 
that the high-energy emission originate from the same internal 
dissipation process as that of the sub-MeV emission (Abdo et 
al. 2009a). Taken together, these results suggest that high- 
energy emission of GRBs detected by Fermi/LAT during the 
prompt bursting phase is dominated by the high-energy emis- 
sion arising from the internal origin, rather than the onset of 
the external shock, although the external forward shock can 
readily account for the extended high-energy after the prompt 
phase (Kumar & Barniol Duran 2009, 2010; Ghirlanda et al. 
2010). 
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APPENDIX 

Appendix: The forward-reverse shock scenario in different spectral regimes 

As we mentioned in §2, there are three different regimes in which the light curve in the radio band is affected by the synchrotron 
self-absorption effect and thus its behavior is different from that in the optical band at 4.67 x 10 5 s. In fact, these three cases can 
be categorized into two cases, i.e. (1) v ar > v m and v M > ^ ra( jio, (2) v w < and v m > ^ ra dio- We take £54 = 2.0 and n_3 = 0.40, 
as constrained by the late radio, optical and X-ray emission and a typical value of p = 2.5 for the reverse shock is used here. 

THICK SHELL CASE 

The necessary condition for the thick shell is 790 > 7d ec , so we get 773 > 1.91. According to equations (3), (19) and (20), we 
have 

V W = 1.05 X lOV^l^r 2 -!"-^ 2 ( A1 ) 

a = 5.13 x 10 7 7/3 6 £ 5 4n:3 /2 e; r 5 -ieBy-i ( A2 ) 



and 



4 Using a numerical code, He et al. (2010) find that a single jet model can j ata b e f ore ^ j s 

explain the broadband afterglow of GRB090510, except the early time LAT 



at 4.67 x 10 5 s. So the SSA frequency is 

/ 4.52 x 10 V /5 4Vf e^^Hz^ < 1 
V " = \ 2.52 x 10«, 3 2/13 < /13 „!f 4£4£Hz >fl > 1 ( } 

far > Vmr and V w > V mdio 

According to the second equation of (A3), we get vl^^^-^Bi 2 -! > 37.4 or 773 > 1 .67 x lO 10 ^^^^ by requiring u m > v m d{ . 
Such a large value of the initial Lorentz factor is unlikely, so we can exclude this case. 

far < v mr and v M > v mdio 

By requiring v w < v ml , we obtain vt^sl^-i^-et -i e B^-i > 5.13 x 10 7 . On the other hand, from v- dI > v m dio we obtain 

?73^ 5 £ , 5 4^ 5 «l3^ 5 e er _ieg{ 5 _ 1 < 5.32. Combining these two inequalities together and taking £54 = 2.0 and n_3 = 0.40, we have 

773 < 9.2 x 10~ 3 eg{ 4 _j. Such a low initial Lorentz factor is inconsistent with 773 > 1.91, derived from the pre-condition for the 
thick shell case, so this case can be also excluded. 

THIN SHELL CASE 

The thin shell case, on the contrary, requires 7d ec > 25s, or 773 < 1.91. On the other hand, the fact that prompt high-energy 
photons can escape from the source without annihilation with low-energy photons puts a low limit on the initial bulk Lorentz 
factor, which should be larger than a few hundreds. Another constraint comes from the high-energy afterglow emission. If we 
attribute the high-energy emission after 170 s detected by LAT to the forward shock, the deceleration time should be shorter than 
170s and then we get 773 > 0.92. In the thin shell case, we have 

Vw = 6.43 x 10V4 /2 4-i e Br 2 -i Hz , (Bl) 

a « 6.72 x l^'%l l7A nT^U4J-, (B2) 



and 



at 4.67 x 10 5 s respectively. 



2.02 x 10 V 5 4M 5 <-i4£!Hz,a < 1 

3.99 x 10V V13 4 /13 «-f e^e^Hz,,, > 1 ^ 



T^ar ^ ^mr ? far ^ f radio 

From^ ar < v ml , we get 6.72 x ^r/^^EJ^n^^e^^e^/^ < 1, while from v w > v m&0 , we get 773^ 'E^n 2 ^ 'eer,-i4{f-i > °- 42 - 
Combining these two inequalities together, we have 77 3 8/3 < 1.46 x lO^f-i or m < 0.036e B { 1 ^. Since such a low initial bulk 
Lorentz factor is inconsistent with the low limit on the initial Lorentz factor constrained by the high-energy gamma-ray photons, 
this case can be excluded. 



jy- dv > 7^ mr , V dx > f radio 

Condition v w > v m dio requires ri^^^E^^r^J^e^^e^^ > 21-3. If we substitute £54 = 2.0 and n_3 = 0.40 into this inequality, 

we get % U ^ 3 ^-i^-i > 19-1 or 773 < Q.Q2>\e 6 J r l \e^\. Even with e er -i = 10 and e Br -i = 10, we still have 773 < 0.28. Due to the 
same reason as the previous case, this case can be ruled out either. 

Consequently, we conclude that either v m or v mx cannot be above the radio frequency at the second radio detection time 
4.67 x 10 5 . The light curve in the radio band at 4.67 x 10 5 s should behavior the same way as the optical emission around 5000s. 
And based on our discussion in §2, we conclude that the FS-RS model is not viable for the broadband data of GRB 090902B. 
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Table 1 

Forward-shock parameters in two-component jet model. 
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Table 2 

Reverse-shock parameters in two-component jet model. 
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